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Abstract 

We derive the analogue of the QCD low energy theorems for the scalar and pseu- 
doscalar gluonic correlators in nuclear matter. We find that the scalar correlations 
are depleted while the pseudoscalar correlations are enhanced to leading order in 
the nuclear matter density. We briefly discuss the consequences of these findings on 
the QCD spectrum. 
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1 INTRODUCTION 

QCD Low energy theorems (LET) for the gluonic correlation functions have 
provided insights into the QCD dynamics both in the vacuum and in the 
nuclear medium. Indeed, the LET for the pseudoscalar gluon correlation func- 
tion provides a powerful constraint on the contribution of the rf mass [1,2] 
and the vacuum gluon condensate [3]. We recall that a quantitative calcula- 
tion of glueball masses using QCD sum rule analysis [4, 5], require the use of 
LET to provide the necessary subtraction constants needed in the dispersion 
analysis [5, 6]. Recently, the extension of the LET to finite temperature and 
density have been derived for the correlation functions of the scalar gluon op- 
erators^]. These theorems provide strong constraints to be satisfied by lattice 
calculations or any effective model calculations of QCD in matter. Here, we 
will derive anew the LET for the scalar gluon correlation functions at low 
density and show the equivalence of our result to that of ref. [7]. Then, we will 
derive similar LET for the pseudoscalar gluonic correlation functions in two 
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ways : First, by assuming that the self-dual gauge configurations are domi- 
nant in vacuum and (dilute) matter; Second, by using a generic heavy-quark 
expansion. In both cases, we will derive general relations on the density depen- 
dence of the correlation function in terms of differential operators. We briefly 
comment on the relevance of our results to the QCD spectrum. 



2 LET in vacuum: review 



Let us start with the definition of the scalar and pseudo scalar gluonic two 
point function. 



S(Q 2 )=t [ d 4 xe^{T*^G 2 (x)^G 2 (0)) 

J 47T 47T 

P(Q 2 )=i J d 4 xe iqx (T*^GG(x)^GG(0)}. 



Here, G 2 = G^G^ and GG = G^^e^apG^. The low energy theorem for 
the scalar gluon operator follows from the formula[3], 



d 



(O)=tJd'x(T*O(x)g 2 G 2 (0)), (2) 



where go is the bare coupling constant of QCD. Now using renormalization 
group arguments, one notes that the bare coupling is related to the ultraviolet 
cut-off Mq via, 



(0) = const 



M exp(- 



Sn 21 " 



(3) 
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with 6=11 — |jV/ to one-loop. Therefore, the left hand side of eq.(2) yields, 



(0)=d™£(0) (4) 



Substituting O = ^G 2 (0) into eq.(2) and using eq.(4) we find the LET for 
the scalar gluon correlation function. 



S(Q 2 = 0) = ^G 2 } (5) 
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The situation is more subtle for the pseudoscalar gluonic correlation function. 
Indeed, in the presence of light quarks, the pseudoscalar gluon field is a total 
differential of the light quark axial- current, i.e. ^-GG = S 9 ^?7/*7 5 ? ( m q ~* 
0). Hence, 



P{Q 2 = 0) = . 



(6) 



This is not true in the absence of light quarks (NLQ). The approximate LET 
for this case was originally derived in ref.[3] by assuming that self-dual field 
configurations dominate the functional integral [3] over the gauge fields. To 
show this, consider the QCD functional integral with nonzero 9, 



Z = j [dA]exp(U) 



(7) 



where 



d x 



4 32vr 2 



(8) 



Differentiating eq.(7) twice with respect to 9, we have, 



-i J d 4 x(T* 



2 2 

% GG(x),^GG(0)) 



327T 2 



327T 2 



(9) 



where e is the vacuum energy density. To evaluate the left hand side of eq.(9), 
we note that the energy density is related to the trace of the energy momentum 
tensor 9^, 



4^ ^ v 16a s 1 x 32tt 1 



(10) 



Now we assume that the self-dual (anti-self dual) field dominates the functional 
integral. This changes the euclidean action as follows, 



IE 



I 



d 4 x 



32tt 2 



d 4 x 



V + 10 )G> 



'ID 



where G is obtained by redefining the gauge field — > j^A^- Eq.(ll) implies 
that within this approximation, the ultraviolet cutoff M is related to 9 in a 
physical quantity by a simple replacement of eq.(3), 



3 



d 

(12) 

for small 9 1 . With eq.(12) and eq.(10), it is straightforward to calculate 
the left hand side of eq.(9). We obtain the LET for the pseudoscalar gluonic 
correlations [3], 

Pnl Q (Q 2 = 0) = S(Q 2 = 0) = -y (^G 2 ), (13) 

where the superscript sd refers to the fact that it is based on the self-dual 
assumption detailed above. 

We recall that the difference between the case with light quarks and without 
light quarks is the contribution due to the 7/ [1-3], which adds up to give the 
required zero, i.e. 



(O) = const 



P(Q 2 = 0) = P NLQ (Q 2 = 0) + ^ , lV)l = 0. (14) 

One should note that similar LET for higher point gluonic correlations can be 
obtained by taking further derivatives either with respect to 9 or —1/4^0- 



3 LET for the pseudoscalar gluonic current: Heavy Quark Mass 
expansion 



The above derivation for the pseudo scalar current was obtained within the 
self-dual approximation in the path integral. This approximation can be averted 
Indeed, let us consider a regularization scheme with 9 7^ in the presence of a 
heavy quark with mass toq. Substituting the pseudoscalar gluon field for the 
operator O in eq.(2), we have, 

-^(^GG) = i J d 4 x{T* ^ G G W ^(0)> . (15, 
Now, we will make use of the heavy quark mass expansion[8], 

m Q Qi'y 5 Q = ~GG + o(^ T ) 
8tt \m 2 Q J 

1 We note that for large 6, 2ir periodicity is recovered through the branch-structure 
of e. 
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Substituting eq.(16) into eq.(15) we obtain, 

d(-l/4 g g) (6mq< ^ 75g) = 1 1 ^ G ^H87T 2 m Q QQ(0) ) . (17) 

In general, such substitution is not allowed when the right hand side of eq.(17) 
has an e tqx factor as one is calculating the operator product expansion. How- 
ever, since we are looking at the LET dominated by non-perturbative space- 
like field configuration [9], such substitution is valid. 

Consider now making a chiral rotation in 9 — > 9 + f . This will only affect the 
heavy quark mass by mq — ■> mge* 7571 "/ 2 . As a result, the scalar quark condensate 
turns to the pseudoscalar quark condensate and vice versa. Substituting the 
heavy quark operators back to the gluon operators, we obtain the following, 



'— G 2 ) = i [ d 4 x(T* ^GG(x)^GG(0) > . (18) 
ir J Air Air 



Making use of eq.(4), yields 



P h N 9 LQ (Q 2 = 0) = ~{f G 2 ) , (19) 

where, the superscript hq refers to the heavy quark mass expansion used. This 
is different from eq.(13) as obtained within the self-dual approximation. Given 
the uncertainty in the value of the non perturbative gluon condensate, both 
results are consistent with the vacuum phenomenology in eq.(14). 



4 LET for the scalar gluonic current in nuclear matter 



Here we derive the LET for the scalar gluonic current at low density. The 
same derivation holds for finite temperature. The starting point is taking the 
expectation value of eq.(2) at finite density. Then we make a low density 
expansion of the left hand side, 

SFlkl) {0) > = dTTM) + + °<'» ■ < 20 > 
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Here, {-) p denotes the expectation value at finite density, p is the nucleon den- 
sity and \N) a nucleon state normalized as (N\\N) = (27i) 3 E^/m]\fS 3 (0). This 
is the linear density approximation. The derivative of the nucleon expectation 
value is simply obtained by substituting d — > d — 3 in eq.(4). The derivative of 
the nucleon density can be obtained from the fact that the chemical potential 
p = p 2 f/2rriN is an external parameter and is independent of go- That is, using 



d ^'Uo (21) 



<i(-l/4 Jo 2 ) \2m„ 



and 



we find, 



d 32tt 2 

m N = ——m N , (22) 



d o 327r 



2 



di-l/HY* b 



P f = —f—2m N p, (23) 



d 32tt 2 3 . . 
2TP = P^-o- (24) 



d(-l/4giy ^ b 2 
Therefore, 

d 32tt 2 / _ 3 



rf(-l/^ 2 ) 



(0) p = (d(O) + (d- l)p(N\0\N) + 0{p)) . (25) 



From this, it follows that the changes of the correlator for the scalar gluonic 
current to leading order in density is, 



AS(Q 2 ) = S(Q 2 ; p) - S(Q 2 ; p = 0) = ^-p(N\^G 2 \N), (26) 

where S(Q 2 ,p) is the correlation function defined in eq.(l) calculated at finite 
density p. 

We note that eq.(25) is also consistent with the result of ref. [7] applied to 
leading order in density. In ref. [7] it was noted that a physical quantity of 
mass scale d has the following dependence on the temperature or density. 



(0) P ,T = A d f(j,^), (27) 
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where (-) p ,t is the expectation value at finite density and temperature, and the 
mass scale A depends on the coupling as in eq.(3) with d—1. Then, it follows 
that the differential operator acting on physical parameter can be replaced as, 



' -(%=?^4-4l(% (28) 



d(-l/4gi) x /p ' 2 b \ dT HJ d/i y 

When we apply this differential operator to our eq.(20) and note that /ij^ 
|p, we recover eq.(25). 



5 LET for the pseudo-scalar gluonic current in nuclear medium: 
self-dual approximation 



The LET for the pseudo-scalar gluonic currents at low density can be similarly 
obtained in the self-dual approximation. Here, we consider the case with no sea 
light quarks, stressing the effects of matter on the gauge configurations. The 
starting point is the generalization of eq.(9) to finite density. The difference 
compared to the scalar case is in taking the derivative with respect to 9 instead 
of <7q. Again, using the linear density approximation, we have 



^(°)p = ^2 «°>o + P(N\0\N) + O (p)) . (29) 

To evaluate the derivatives with respect to 9, we only need the relation in 
eq.(12) and the fact that dp/d9 = 0. To leading order in density, we have 



dp 3 i 

Therefore, we have 

^- 2 (0) P = (y f{d 2 (O) + ( 9 - + 3(d - 3) + (d - 3) 2 ) P (N\0\N) 

+0(pj). (31) 
Substituting this into eq.(9) gives, 

^Pnlq(Q 2 = 0) = P s N d LQ {Q 2 = 0; p) - P s N d LQ {Q 2 = 0; 0) 
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(32) 



We can also derive a general formula for eq.(31) in terms of differential opera- 
tors involving the derivatives with respect to the temperature and/or chemical 
potential. Again, we have the 9 dependence of a physical parameter through 
eq.(27) and also eq.(12). Therefore, we have 

^(0),T-(iy(l-T^-^J(0) e „ (33) 

which is consistent with eq.(31). In the linear density approximation eq.(33) 
can be written as, 

^«0>o + P (N\0\N)) = (V(O)o + (d - 3 -Tp(N\0\N))) . (34) 



6 LET for the pseudo-scalar gluonic current in nuclear medium: 
Heavy quark expansion 

The generalization of the low energy theorem for the pseudo-scalar current 
based on the heavy quark expansion can be obtained in the same way as 
in the scalar case. The derivation provides for a check on the self-duality 
assumption. Starting from the expectation value of eq.(18) at finite density, 
and substituting eq.(25) into the left hand side, we obtain, 

AP£l (Q 2 = 0) = P h N %{Q* = 0; p) - P h N \ Q {Q* = 0; 0) 

= - 5 -p(N\^G 2 \N). (35) 
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Note that the overall factor is different from eq.(32), but the sign is the same. 
The relation in terms of differential operators can be obtained in the same way 
as in the scalar case by substituting eq.(28) into the left hand side of eq.(18). 



7 Summary 

We have derived the LET for scalar and pseudoscalar gluonic correlation func- 
tions at low density. Although, we have restricted our discussions to two point 
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functions, we can easily generalize this method to higher point functions and 
to finite temperature. Since eq.(14) is still valid in the medium, we expect, 



A 




) 



&Pnlq(Q 2 = 0) < . 



(36) 



The magnitude of the right hand side depends on the approximation scheme 
we have used (self-duality or heavy quark). However, in both cases, the sign 
is negative. To leading order in the matter density and if we were to assume 
that residues do not change, then the rj' mass will increase, while the scalar 
glueball mass will decrease. These observations are worth testing using QCD 
lattice simulations. The faith of the rj' mass in matter may have interesting 
consequences on low mass dilepton emissions, pion thresholds as well as the 
strong CP problem in matter. 
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